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Abstract. Starting from the Oppenheimer-Snyder solution for gravitational 
collapse, we show by putting it into the harmonic coordinates, for which the dis- 
tant Riemann metric is galilean, that the final state of collapse for a collapsed 
star of any mass, including the one thought to occupy the centre of our galaxy, 
has a finite radius roughly equal to its Schwarzschild radius. By applying an 
expression for the gravitational energy tensor, we are able to explain the con- 
centration of stellar material in a thin shell close to the surface, which gives an 
explanation for why such a star does not undergo further collapse to a black 
hole. The interior of the star is characterized by a low density of the original 
stellar material, but, far from being empty, this region is occupied by a very 
high density of gravitational energy; this density is negative and the consequent 
repulsion is what produces the surface concentration of stellar material. 

1 Introduction 



There is a widespread belief that space tells matter how to move [I]. This, we 
believe, has resulted in a profound misunderstanding of gravity. Although the 
error can be traced back all the way to the founder of the modern theory [5], 
there is ample evidence that he nevertheless had a strong inclination to go in 
the direction which we are advocating here. 

Throughout the decade in which Einstein discovered what he called General 
Relativity (GR), he repeatedly attempted to cast gravity within the context of 
a classical field theory by constructing a tensor for the energy and momentum 
carried by its field[3]. Perhaps the culmination of this effort was the article[4] 
in which he derived the quadrupole formula for the emission of gravitational 
radiation from a bounded nonspherical source. He had already changed his 
mind twice about such radiation when he wrote that article; throughout the 
subsequent six decades, not only would he himself change two more times, but 
his oscillations of opinion would be refiected in the communal understanding 
of gravity, causing confusion among virtually all of the leading scholars in the 
field |S]. The communal view changed radically when observation of the Hulse- 
Taylor pulsar [BJ confirmed that this system is losing energy at the rate predicted 
by Einstein's quadrupole formula 

But did it? In our opinion the change is incomplete and inadequate. Proba- 
bly most of us now confidently expect that the LIGO experiment [7] will reveal 
that gravitational waves passing through space cause a change in the path of 
light signals exchanged across that same space. The light signals are items of 
"matter", which is being moved but by what? We have become lazy in repeat- 
ing, as a kind of mantra, that the moving agent is space itself. That may be an 
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adequate description of the static solar field's effect on planets and passing light 
signals, but the time has now come to distinguish between the gravitational 
field and the space carrying it. The field has energy which may be localized in 
a determinate manner; the energy has mass which itself gravitates, and above 
all it must be described by an energy tensor. 

In spite of appearing in Einstein's own articles, both before and after GR was 
conceived [3], the energy " pseudotensor" has not achieved full tensorial status 
right up to the present day. There is a big obstacle to be overcome, and it was 
pointed out by Hilbert[H], who was the codiscoverer of the basic Hilbert-Einstein 
field equation of GR. This is that the metric and its curvature are the only proper 
tensors which may be constructed out of 17^1/. The analogy with classical field 
theories suggests that an energy tensor be formed from the derivatives of g^^i,, 
but normal derivation is not tensorial and covariant derivation in the Riemann 
metric gives zero. 

It must have been physical intuition which caused Einstein to ignore Hilbert's 
objection and continue with his "pseudotensor" to deduce the quadrupole for- 
mula. In order to make his argument work he made a particular coordinate 
choice, using what is now called the harmonic system, thereby violating the 
Principle of Equivalence (PE) , which he had put at the centre of his derivation 
of the Hilbert-Einstein equation, and which states that no particular coordinate 
system is to be preferred. Now we are able to see that, rather than introducing 
a favoured frame, what he was actually doing was to recognize that the gravi- 
tational field is carried by the familiar Minkowski space of so called "Special" 
relativity. The notion of the gravitational field as being like the electromagnetic 
field of Faraday and Maxwell was taken up by Einstein's contemporaries de Bon- 
der and Lanczos[Tn][ll], and subsequently developed by Fock[l2j, Rosen[T3j and 
Weinberg[T3]. This field interpretation of gravity has subsequently been devel- 
oped, by Logunov and Mestvirishvili[T5], and by Babak and Grishchuk[T6j, to 
the point where the Minkowski metric is explicitly present in the field equations. 

In advocating the necessity for including the Minkowski metric in the field 
equations, we will also be stressing the need to distinguish between the various 
requirements which have been laid on the theories of gravitation, namely co- 
variance, gauge invariance and the Principle of Equivalence. We shall argue for 
maintaining the first, abandoning of the second, and accepting only the weakest 
form of the third, in the form of the Eotvos Principle (EoP). 

Our new contribution to the field theory of gravity is in gravitational collapse. 
Starting from the Oppenheimer-Snyder solution[17j. we show by putting it into 
the harmonic coordinates, for which the distant Riemann metric is galilean, 
that the final state of collapse, for a star of any mass, including the one thought 
to occupy the centre of our galaxy, has a finite radius roughly equal to its 
Schwarzschild radius. 
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2 The energy tensor 

We follow Babak and Grislichuk[l6] (BG) with some changes of notation. The 
self-interacting tensor field h"^ {x) is defined on an underlying (fiat) Minkowski 
space with metric 

d(7^ = 7^j, (x) dx^dxi, . (1) 
The d'Alembertian operator is 

a^^^^-'D^D, , (2) 

where 

r''7.x = 6^^ , (3) 

and Df^ denotes covariant differentiation with respect to x^ in the Minkowski 
metric. We define 

and we shall denote covariant differentiation of $ simply by a lower index, that 
is 

= D^<i>"^ ^ 7f/ + hff = hfj' . (5) 
We define also the inverse field tensor 'i'ajs by 

^a0'i>^''=Sl . (6) 

The field equations are determined from a lagrangian density 

L = Lg + L^ . (7) 

The gravitational lagrangian is 
where 

(9) 

and 

$ = det ($°^) . (10) 

This is equivalent to an expression first given by Fock (Ref[T2], Appendix B). 
The "matter" lagrangian L„i is a function of all the nongravitational (particle 
and field) quantities, collectively labelled by and of h"^ , which enters only 
through the combination ^—7$"'^, that is 

im-im \/=7*"'') ■ (11) 
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The field equations are then obtained by varying L with respect to /i"'^, or 
equivalently and are 

2k, 5Lg_ _ _ „ pa _ pa pl3 , 1 pa p/J _ „^ 



BG showed that these are equivalent to the standard Hilbert-Einstein (HE) 
equations of GR. First one makes the identification of the inverse field tensor 
with the Riemannian metric tensor, namely 

^?a/3 = ^^'a/3, 7 = det(7"'5) , (13) 
from which it follows that 

det(g„^) =.9 = $/72 . (14) 
Then, in terms of the new variables, one finds that 



the right side being the contracted curvature tensor of the Riemannian metric. 
On the other hand, defining the inverse metric tensor G^" (conventional notation 
g^") by 

G'^Va^^^ , (16) 
and the material stress tensor by 



we obtain, since G^^" = ^7/$$^'', 



= T^iu ~ ■7:g^ivG°'^Tafi = T^^ - -g^ivT . (18) 



With these substitutions (fT2|) becomes 

1 



^Mi- = «c 1^7),,. - -gf_iuTj , (19) 

which is the Hilbert-Einstein equatiorQ. 

The energy balance is obtained by varying L with respect to 7^^, that is 
(note that this variation takes account of 7'"' occurring in $^''), 

, dL^ _^ ^20) 



^Our derivation follows a long line of field theoretic ones, starting with Hilbert, who used 
the lagrangian R^—g. The lagrangian of BG is a covariant version of one used from the years 
shortly after the establishing of GR (see, for example, Eddington^ section 58). 
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The first term is given by 



(21) 

where 



an.' 



Hence, making use of the field equation ^I2\i , this first term becomes 



The second term is given by 



a 

(24) 



^^"^ Vr'^-VHS^ ■ (25) 



Hence (1201) becomes 



2k " " ^ ' s ^ y 2 J d^^^P 

= + ^ (^Q<-^^QP- _ l^^^G'^''^ (^T,,, - ig^.G'^-T,,^ . (26) 

In the usual notation the inverse Riemann metric is G"'' = g"^, and it is the 
tensor used for raising indices, so that the right side simplifies to give 

Da^Dp ($''''$"'3 _ ^M"^''/^) ^ 2Kt^"', t^"" = {t^'/ + • (27) 

Now the expression DaDpD^ {^t^^i^^P _ jg antisymmetric with re- 

spect to a /ia interchange, from which we may deduce the covariant conservation 
equation 

Dpt^"" = . (28) 

The theory summarized in this section has all three basic properties mentioned 
in the previous section, that is covariance, gauge invariance and weak equiva- 
lence as summed up in EoP. The latter may be briefly stated as "All forms of 



^This expression was given in noncovariant form, that is with ordinary instead of covariant 
derivatives, by Landau and Lifshitz|18| and also by Fock|12|. as the gravitational energy 
"pseudotensor" . 
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energy gravitate equally". It has an active aspect characterized by the single 
coupling constant k in (j27l) . and also a passive aspect, contained in the equation 

y^j^f,^ ^ . (29) 

Of course, the latter property is a consequence of together with the Bianchi 
identity 

(r'^" - ^g'^'R] - . (30) 



3 The harmonic condition 

The condition 

= (31) 

has a history almost as long as GR itself. It was used by Einstein[4] in his 
derivation of the quadrupole formula, but was almost immediately criticized, 
for example by Eddington (Ref[5, page 130), for being noncovariant; Edding- 
ton said that gravitational waves "travel at the speed of thought", and his 
criticism resounded throughout the subsequent six decades, finding fame in the 
title of a seminal history[5 of gravitational waves. In the first decade of GR the 
harmonic condition was developed by de Ponder jlOj. and also by Einstein's re- 
search assistant Lanczos[TT|: later it was championed by Fock[12 , who claimed 
that this choice of coordinates was necessary in order to guarantee the condi- 
tion that there are no ingoing waves at large distance from an isolated massive 
object. Fock, and also Weinberg[14] have shown that this condition, as well as 
being an essential part of the field interpretation of GR, is extremely practical 
in all analyses of the far fields of such a system, for example in establishing 
post-Newtonian expansions. Subsequently [TH] it was argued by Logunov and 
Mestvirishvili that, in its covariant form 

^ , (32) 

the harmonic condition, far from being just a choice of coordinate system or 
"gauge" , is, along with the Hilbert-Einstein equation dealt with in the previous 
section, an essential field equation of gravitation. These authors showed that, 
by including the Minkowski metric "f^^, explicitly in a lagrangian 

Lgi = A (2V^7^,a>^'' - V^) , (33) 

so that the total lagrangian is 

L = Lg+Lgl+ L„, , (34) 

the covariant harmonic condition does indeed appear as a field equation. Such 
a modification causes the Hilbert-Einstein equation to acquire two additional 
terms, one of which is the familiar " cosmological constant" term (what Einstein 
termed his "biggest mistake", which, however, has now come back into fashion). 
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The cosmological effect of tliese two terms has been explored [T^] with A related 
to the Hubble constant, but in this article we shall assume that A is negligi- 
ble, so that the Hilbert-Einstein equation remains unmodified. The extra term 
5Lgi/5h^^^ in the field equation violates the gauge invariance, though not, we 
stress, the covariance. In general this additional term may result also in the 
violation of EoP, in the form (P^ . but the imposition of the field equation (|32p 
is precisely what prevents this occurrence. It should be noted that the loss of 
gauge invariance arises precisely because of this additional field equation, and 
that in field theoretic terms it is more accurate to state that what has been lost 
is the gauge ambiguity of GR. 

As far as the energy tensor is concerned, we may incorporate the new equa- 
tion and write the energy balance (|27l) as 

t^'' = t^^ + i-T^'- , (35) 

where 

+A<^>"'^^p^'^/<^Pp'' . (36) 

Although the two terms on the right side of (1551) may be referred to loosely as 
the "field" and "matter" parts of the energy tensor, it should be remembered 
that they are not just the functional derivatives of Lg and Lm, as given by (j24p 
and (I25|) : some terms were transferred from one of these to the other. For most 
purposes it will suffice to consider them both together, and then they may be 
computed as the left side of ([35|) . that is 

t>'j^ = —DaDfs ($'^''$"'3 _ $M"$^/3) _ (37) 

At this point we specialize to the case of cartesian harmonic coordinates, to 
be discussed in more detail in the next Section, and for which the Minkowski 
covariant derivatives are ordinary partial ones. Then the 00-component of t'^ 
reduces to 

40 ^ J_Q^Q^ ($00^«j _ ^Oz^Oj) ^ (3g) 

which may be expressed as a 3-divergence 

4° = divt, U = e>^J = dj ($0"$*^ - $0'$"J) . (39) 

ZiK 

A general spherosymmetric field may be written, in cartesian coordinates 

xq — t. x\ + xi^ — r,Xi/r — Tii, as 

$™ = $1 (r, t) , $0' = $2 (r, t) n„ = $3 (r, t) n,n, + $4 (r, t) (% - n,n, ) 

(40) 
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and its determinant is 



$ = ($1$3 - $2) *4 , 



and so the cartesian tensor on the right side of becomes 



that is 



t = 



$1$4 (^,, 



(41) 
(42) 
(43) 

1-4 ' \-*'4 /J 

As a first apphcation we use to find the field energy distribution in the 
exterior of a nonstatic spherosymmetric T^^'^ with total mass M . By Birkhoff 's 
theorem the field is static, and is the harmonic version of the Schwarzschild 
metric, that is (see Ref[I2j 209-215 ) putting m ^ GM/c^ 



2k 



' <i> 2 / $ 



$1 = 



(r + m) 



(r — m) ' 
it follows that 



$2 =0, $3 = 



$4 = -1, $ = - 



Mc^n / r + m,\^ 2r — m 
Airr'^ \ r J 2r — 2m 



We may define a new vector ti by subtracting a divergence-free vector 

Mc^n 



ti = t - 



(44) 
(45) 

(46) 



and then the total (matter plus field) energy contained in r > ri is Mc^^ (^i)j 
where 



divtid-^r = 1 



r>ri 



ri + m\ 2ri — m 



2ri — 2m 



(47) 



but in this case it is entirely field energy. Note that this quantity is negative; if 
ri ^ m it is 



— + o( — 



(48) 



and, for example, the absolute value of the field energy in the exterior space 
of our Sun, for which m/ri « 2 x 10^^, is about three times the mass of our 
planet. 

More generally, we now show that, in the spherosymmetric case, there is a 
simple connection between the cartesian tensor Qij and the spatial components 
Qij, of the Riemannian metric of GR. In terms of the components of the field 
tensor we obtained (see equation (221)) 



(49) 
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The inverse field tensor is 



00 



3^4 



$ + ^ (% - niUj) (50) 



Since in these coordinates 7 = — 1, the corresponding Riemannian metric tensor, 
from p^ . is 



gap = \/-**a/3 , 



that is 



500 



9oi = 



9^J 



$1$ 



Converting to spherical polar coordinates we have the metric 



ds^ = goodt'' + 2gordtdr + grrdr^ + gee {dO^ + sin^ 



where 



90r 



='2'^4 



$4 



Then (1421 becomes 



which gives 



= divt, t 



n / 9.gee .gee - gr 



K \ dr 



gee 



(51) 



riin,- H — (fc - rijn,) 



(52) 

(53) 

(54) 
(55) 
(56) 



and hence the total energy (matter plus field) in r > ri at time Hs /i (ri , t) Mc^ , 
where 

gee 



2mr^ 



, 2 

+ r grr - r 



dr 



(57) 



4 The Oppenheimer-Snyder solution 



We now extend the latter calculation to a continuous mass distribution in order 
to discuss the gravitational collapse of a spherical star. In contrast with the 
negative energy distribution of the field, the mass tensoiH T^'' has the positive- 
energy property that its contraction T^'^'q^rj^, with any timelike covariant vector 
T}^ is positive. One of the few exact solutions for a Riemannian metric derived 

^We use the short name, favoured by Fock, for this tensor. Neither it nor the more usual 
"material stress tensor" are quite exact according to our present perspective, because the 
gravitational field is material and has both energy and mass. 
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from such a mass tensor is that of Oppenheimer and Snyder[T7] (OS) derived 
in a comoving and synchronous coordinate system, for which the only nonzero 
component of T'^'^ is a positive T"", 



ds^ = dr' ~ {d*jiWf dR^ - W 



^2 /-J^s ^sin^^-^•^2^ 



(58) 



where the operator d'^ denotes differentiation with respect to R at constant t, 
and 

IW 3r^/2^ 



W = u 



(59) 



The function u (R) gives the steUar material density, and is a positive monotone 
increasing function for i? > with u (0) = and u(oo) — 1. We now choose 
units such that 2m = 1 and define a new coordinate 



IR^ 3t 



1/3 



to replace r, so that, denoting differentiation at constant v by Or 



3/2/1 / 

'r^ ~^ 



and 



Then the Riemannian interval in terms of (_R, u, 6, 0) is 



ds^ = {PdR - 2v^dvf - ev^dR' - ^-^^2 , ..„2 ^..2^ 



+ sin^ 



that is 



CSC 



2 4 
-U V 



(60) 

(61) 
(62) 

(63) 
(64) 



The contravariant inverse of this (we are now reverting to conventional notation) 
is, 



-,RR 
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0" 



r,VB. 



1 



2Cv 



2„6' 



2„4 



(65) 



The Cartesian harmonic coordinates are obtained, (see Refjl4l ppl65-168), 
as solutions of the equations 



□ga;^ = (/i = 0,1,2,3) , 



(66) 
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where Og is the Riemannian (not to be confused with the Minkowskian ^) 
d'Alembertian 

□g = -La, (g'^^V^dp) . (67) 



In spherical coordinates {t = xo,r = ^yx1 + X2 + x'^) these reduce to two equa- 
tions for what OS term the "exterior time" and "exterior radius", namely 

Qt = 0, (Q-^)^ = ' (68) 
where Q — —£,'^v^Dg and is given by 



9R + i,d?l^[dn + ^,d.)~^,d.v'^d.. . (69) 



However, in the present treatment the coordinates (t, r) describe both the exte- 
rior and the interior of the star. 

The free-space limit of the last Section corresponds to u = 1,(3 — VR,£, = 
y/R/v^, for which the solutions of the latter equations are 15 

r = «2_ t^-R^I^ -2^{v), ^{v)^\ (w^>\) , (70) 

^ -J JwQ ^ ^ 

and, substituted in the OS metric ((63)) . these give the harmonic version of the 
Schwarzschild metric, that is 

2_2r-l 2 2r + l 2 ( „ , ^\ fjQ2 , „• 2 /, ,,2 



ds' = ^^dt' - - (^r + - j {de'+ sin' edcp') , (71) 

and this corresponds to the far field ^'"^ we used in the previous Section. Both 
of the PDEs ((68|) have a singularity along the curve v = vq (R), specified by 



dvo 1 , ^ 
"0 



In our earlier article pO]. we showed, for a particular choice of u, how to solve 
(|68p numerically by integrating along a family of characteristics, satisfying the 
ordinary differential equation 

The operation of differentiation along such a characteristic is 

d^^dn-(^-^)d. , (74) 



2 2v^ 

and the operator Q may be expressed as 
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where 



2w' 



5^ 
2^ 



(76) 



We now define the physical region as R > l,v > vq (R); it is our contention 
that ((68)l describes the whole evolution of the collapsing system, without any 
singularity, in this region, and the final stage of the collapse, corresponding to 
the limit t {R, v) — > +oo, is described by the values of r and t close to the 
limit characteristic v — vq (R). A necessary preliminary to a numerical study 
of these equations is the obtaining of asymptotic expansions, both for large R 
and for large v when R is small, and this was done in our earlier article [2D]. An 
extension of this calculation may now be made by applying these earlier results 
to the total energy distribution given by (j57p . putting 



and 



dv 



2 4 
-U V 



,dR 



(77) 



\ dr ^ dr 
~J^~dR\dR' 



'dv 



where 

to give 
l-fi{R,v) 



j_d^dr_ d^r dt 
dR dv dR dv 



(78) 



(79) 



^ 2r (2d^t /I 1\ dt 
J \ V dR \u V I dv 



4r2 d^t /d^t dt 

J^u^ dR \dR^^d^ 
(80) 



It is appropriate to mention again here that t^j^ (r, t) , from which the quan- 
tity /i (_R, v) was derived, is indeed a tensor. If, for example, we should choose 
to transform it back to the original OS coordinates, this would be possible, pro- 
vided we bear in mind that in those coordinates the partial derivatives da and 
dp would be replaced by the Minkowski covariant derivatives Da and Dp] the 
latter, of course, are not the same as the Riemann covariant derivatives Vq, and 
V^. In order to describe the actual energy distribution we need /i as a function 
of (r, t), which requires rather complicated interpolations from /i [R, v) , r (i?, v) 
and t {R, v) , and we shall describe such a detailed calculation in another article. 
Nevertheless we venture to draw some qualitative conclusions from the results 
already obtained. 



5 What we expect to find 
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Figure 1: The energy content 1 — /i(r), in units of Mc^, contained within a 
sphere of radius r measured from the centre of a collapsing star, as a function of 
r for fixed t. Between A and B the energy is mainly gravitational and negative, 
between B and C it is mainly stellar material and positive, and beyond C it is 
again gravitational and negative. The radius r is in units of the Schwarzschild 
radius. 

The field part of the 00-componcnt of the energy tensor is negative, as we found 
in an earlier section when we examined the far field. In the deep interior of the 
star we expect to find that this component is substantial, and indeed bigger in 
absolute magnitude than T°°. Being negative, it produces, in a bootstrapping 
manner, its own gravitational field, and this is repulsive. That explains why, 
as we already found in our earlier article [2^, the material energy is con- 
centrated in a shell near the surface as i — >■ +oo. The density function u (R) 
describes a star with an initially diffuse corona, but the formation of a distinct 
surface is an expected part of its evolution. As pointed out above, a substantial 
amount of computation is required in order to obtain a quantitative expres- 
sion for the total (gravitational energy plus stellar material) mass distribution, 
H (r, t) . In Figure 1 we give an artist's impression of this profile, for fixed t, at a 
fairly advanced stage of the collapse process, when the surface shell has begun 
to form. 

The key element in our approach, stressed also in our earlier article [20]. is 
that the coordinate-free topological analysis, now common in GR, fails to ac- 
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knowledge the central role of the energy tensor. We accept the Hilbert-Einstein 
field equations, but add a supplementary set which results in the inertial (that 
is harmonic) coordinate system being privileged. That means maintaining both 
covariance and a weak principle of equivalence, in the form of EoP, but reject- 
ing gauge invariance. As was emphasized first by Fock, this latter " loss" means 
taking the "G" out of GR, leaving us, nevertheless, with a theory which Fock 
himself called Einstein's Theory of Gravitation. 

As far as the currently fashionable notion of "black holes" is concerned, our 
conclusion is that there is a limit to the degree of compactification undergone in 
a collapsing gravitational system, and this limit is set by the properties of the 
gravitational field itself rather than the presence of the other forces of nature, 
like neutron degenerative pressure. Relatively small "condensars" are observed 
as neutron stars, or pulsars, but heavier ones, like the one thought to lie at the 
centre of our galaxy, are more diffuse, probably of white-dwarf density. Our 
analysis is of an idealized, purely gravitational system, and our key finding is 
the limit characteristic curve, which represents the limit t — > +00. This curve 
is equivalent to the "event horizon" popularly acknowledged to come out of 
coordinate-free GR, but for us there is nothing the other side of the limit char- 
acteristic. The other side of our "event horizon" lies outside the physical space; 
to go there is to go "to infinity and beyond" along with maybe Buzz Lightyear 
or Doctor Who! 
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